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Abstract 

A  two-dimensional  thermal  model  was  presented  to  predict  the  temperature  distribution  of  cylindrical  8-Ah  Ni/MH  battery.  Under  the 
forced  convection,  the  temperature  rise  of  the  battery  is  up  to  about  37,  42  and  5 1  °C,  and  the  temperature  profiles  become  non-uniformed 
at  the  end  of  1 C,  2C  and  4C  rate  charge,  respectively.  It  is  indicated  that  the  increase  of  the  convection  coefficient  can  decrease  the  battery 
temperature,  however,  lead  seriously  to  the  less  uniform  temperature  profile  across  the  battery.  The  numerical  studies  indicate  that  the  increase 
of  thermal  conductivity  can  improve  the  uniformity  of  temperature  profile  to  some  extends.  The  battery  temperature  increases  obviously  when 
charged  at  higher  rates.  Overcharge  can  result  in  an  increasingly  higher  temperature  rise  and  a  steeper  temperature  gradient  within  a  battery. 
©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  spirally  wound  design  is  generally  employed  in  the 
manufacture  of  cylindrical  Ni/MH  battery.  This  design  is  of 
importance  to  battery  in  the  improvement  of  the  energy  and 
power  densities  because  of  its  lesser  accessories  when  com¬ 
pared  with  the  prismatic  design  [1],  therefore,  it  is  widely 
adopted  in  the  batteries  for  hybrid  electrical  vehicles  (HEV) 
emphasizing  high  specific  power.  However,  the  specific  area 
is  smaller  for  the  spiral  design,  and  the  cylindrical  batteries 
retain  more  heat  and  less  uninformed  temperature  profiles 
under  the  larger  loads,  seriously  leading  to  the  performance 
fade  in  the  case  of  the  high  power  batteries  for  HEVs.  Bat¬ 
tery  failure  could  also  result  from  excessive  temperature  rise 
and  steeper  temperature  gradient  within  a  battery  [2].  So  it  is 
important  to  study  heat  transport  in  the  batteries. 

Fast  charge  is  technically  needed  in  the  case  of  batter¬ 
ies  for  HEV  applications.  There  exist  lots  of  challenges 
because  of  difficulties  in  controlling  the  charging  protocols 
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under  dynamic  conditions  during  fast  charge.  Thermal  con¬ 
trol  is  one  of  the  challenges  for  fast  charge.  Modeling  and 
simulation  of  batteries  are  useful  methods  for  studying  and 
developing  the  rapid  charging  process  [3].  Thermal  modeling 
based  on  the  first  principles  is  capable  of  providing  valuable 
internal  information  to  help  optimize  the  battery  system  in  a 
cost-effective  manner  [4].  Mathematical  simulation  of  ther¬ 
mal  transport  within  batteries  is  an  effective  tool  to  obtain 
the  knowledge  about  whether  excessive  heat  generated  dur¬ 
ing  the  battery  discharge/charge  process  can  be  removed  and 
how  the  operating  temperature  can  be  controlled  [5].  This  has 
been  done  for  some  cylindrical  batteries — Li/SOCL  battery 
[6]  and  Li/BCX  system  [7].  Based  on  an  assumption  of  uni¬ 
formed  temperature  profile  across  the  whole  battery,  coupled 
thermal/electrochemical  models  were  proposed  to  investi¬ 
gate  a  larger  prismatic  and  a  small  commercial  nickel/metal 
hydride  battery,  respectively  [3,4].  A  two-dimensional  ther¬ 
mal  model  with  a  constant  rate  of  heat  generation  has  been 
presented  to  analyze  a  30  and  a  16  Ah  prismatic  nickel/metal 
hydride  battery  and  packs  [8,9].  In  the  case  of  cylindrical 
Ni/MH  batteries  used  for  HEVs,  however,  there  exist  few 
publications  in  which  present  theoretical  analyses  of  thermal 
investigation.  In  the  present  work,  a  two-dimensional  math- 
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Nomenclature 

A  the  area  for  dissipating  heat  (m2) 

Cbattery  nominal  capacity  (Ah) 

Cp  heat  capacity  (J  kg- 1  K- 1 ) 

Cp,  ref  reference  heat  capacity  (J  kg  1  K  1 ) 
h  convection  heat-transfer  coefficient 

(Win-2  K-1) 

H  battery  height  (mm) 

i  current  (A) 

Zc  thickness  of  the  case  wall  (mm) 

L  total  thickness  of  a  battery  (mm) 

m  as  defined  in  Eq.  (12) 

M  mass  of  battery  (kg) 

q  heat  generation  rate  per  unit  volume 

(J  m“3  s-1) 

Q  heat  generation  rate  (J  s_  1 ) 

r  direction  perpendicular  to  electrodes 

Ro  internal  resistance  ( £2) 

R  battery  radius  (mm) 

t  time  (s) 

T  temperature  of  battery  (°C) 

Too  ambient  temperature  (°C) 

To  initial  temperature  (°C) 

rref  reference  temperature  (°C) 

U  battery  voltage  (V) 

Uo  open-circuit  potential 

V  battery  volume  (m3 ) 

Z  direction  parallel  to  electrodes 

the  volume  fraction  of  positive  electrode 
<5S  the  volume  fraction  of  separator 

<5n  the  volume  fraction  of  negative  electrode 
s  porosity  (%) 

0  state-of-charge 

Xc  thermal  conductivity  of  wall  (W  m  1  K-1) 

An  thermal  conductivity  of  negative  electrode 

(Win”1  K"1) 

Xp  thermal  conductivity  of  positive  electrode 

(Wm-'r1) 

as  thermal  conductivity  of  separator 

(Wm-'r1) 

Xr  thermal  conductivity  in  the  /--direction 

(W  m“ 1  K-1) 

Xz  thermal  conductivity  in  the  ^-direction 

(Win"1  K-1) 

p  average  density  (kg  m-3) 


ematical  model  containing  temporal  and  spatial  coordinates, 
with  variable  heat  generation  rates  and  temperature  gradients 
across  the  battery,  was  presented  for  analyzing  the  thermal 
behavior  and  obtaining  the  internal  temperature  profile  of 
cylindrical  Ni/MH  battery  with  different  design  and  opera¬ 
tional  conditions. 


Fig.  1.  Schematic  diagram  of  a  cylindrical  8- Ah  Ni/MH  battery. 


2.  Model  description 


The  cylindrical  Ni/MH  battery  (as  Fig.  1)  is  composed  of 
a  porous  positive  electrode,  a  separator,  and  a  negative  elec¬ 
trode,  which  were  spirally  wound  and  inserted  into  a  stainless 
steel  case. 

In  general,  a  thermal  model  is  formulated  based  on  the 
thermal  energy  balance  over  a  representative  elementary  vol¬ 
ume  in  a  battery.  For  most  battery  systems,  the  convection 
term  can  generally  be  neglected.  The  following  transient  par¬ 
tial  differential  equation  Eq.  (1)  is  sufficient  to  describe  the 
temperature  profile  in  the  cylindrical  battery  under  the  con¬ 
ditions  of  angular  symmetry  when  the  temperature  gradients 
and  the  heat  fluxes  are  zero  in  the  angular  direction  [4,5] 


3(CPD 

P - AT 


1  3 

r  dr 


+  q 


(1) 


where  T  is  the  temperature,  q  the  volumetric  heat  generation 
rate,  p  represents  the  volume  averaged  density,  Cp  denotes 
the  specific  heat,  Sr  and  Sz  are  the  heat  conductivities  of  r- 
direction  and  ^-direction,  respectively.  For  simplicity,  accord¬ 
ing  to  the  assumption  of  Bernard  et  al.  [10],  the  heat  gener¬ 
ation  rate  is  assumed  to  be  distributed  uniformly  throughout 
the  battery  using  the  volume-averaging  approach  under  the 
assumption  that  no  phase  change  is  involved  and  only  electro¬ 
chemical/intercalation  reactions  exist  for  the  Ni/MH  battery 
of  the  present  investigation,  hence,  heat  generation  rate  can 
be  written  as 


i 


(U  _  u0)  -  T 


dUp 
d  T 


(2) 


where  i  represents  the  current,  V  the  battery  volume,  U 
denotes  the  battery  voltage,  Uq  is  the  open-circuit  potential 
of  a  battery.  The  heat  generation  rate  generally  consists  of 
the  joule  heat  i(U  —  Uq)  (irreversible)  and  the  reversible  heat 
effect  (iT dUo/dT).  Heat  conductivities  are  considered  as  con¬ 
stants  in  the  thermal  model  for  simplicity.  According  to  Ref. 
[1],  the  battery  temperature  coefficient  ( dUo/dT)  of  the  bat- 
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tery  is  presented  as 

d£/o  =  dt^o  _  dC^o 

dr  dr  dr 

where  the  first  term  on  the  right  side  denotes  the  positive 
electrode  temperature  coefficient,  the  second  term  represents 
the  negative  electrode  temperature  coefficient. 

When  a  battery  was  overcharged,  q  consists  of  the  joule 
heat  and  the  oxygen  recombinant  heat  at  the  hydride  elec¬ 
trode,  according  to  Ref.  [11]. 


q  =  1 .48;  +  R^i  (4) 

The  first  term  denotes  the  oxygen  recombinant  heat  at  the 
hydride  electrode,  and  the  second  term  arises  from  ohmic 
heat  from  the  battery  internal  resistance  R^,  which  is  assumed 
to  constantly  be  3  mil  according  to  our  measurement  in  this 
work. 

It  is  further  assumed  that  the  radiation  conduction  from 
the  battery  to  the  ambient  is  neglected  and  only  cooling  by 
convection  is  significant  at  the  boundary,  and  Newton’s  cool¬ 
ing  law  is  applied  to  boundary  conditions.  So  the  boundary 
and  initial  conditions  can  be  written  as  follows, 


dT  (\  /c\_1 

(r-Tio)  at  r=R  (5) 

dT 

— =  0  atr  =  0  (6) 

dr 

dT  (\  lc\~l 

x^=(k  +  xJ  (T-T^  au  =  0  (7) 

~^?  =  G  +  r)  atz  =  H  (») 

T  =  T0  at  t  =  0  (9) 


where  lc  is  the  thickness  of  the  steel  case  wall,  ac  the  heat 
conductivity  of  the  case  wall;  T0 0  the  ambient  temperature, 
and  h  is  the  thermal  transfer  coefficient. 

The  above  boundary  conditions  imply  that  the  core  of  the 
spiral  is  adiabatic,  while  the  outer  boundary  exchanges  heat 
with  the  ambience  through  a  heat-transfer  coefficient  h.  The 
spiral  is  set  at  a  uniform  initial  temperature  Tq  to  which  the 
ambient  temperature  equals. 


Table  1 

Parameters  for  thermal  model  calculation 


Parameters 

Values 

Battery  radius  (mm),  R 

16 

Battery  height  (mm),  H 

60 

Battery  density  (kgm-3),  p 

3900 

Initial  temperature  (°C),  To 

24.4 

Thickness  of  case  wall  (mm),  lc 

5 

Heat  conductivity  of  case  wall  (W  m-1  K-1),  kc 

468 

Positive  electrode  temperature  coefficient  (VK-1), 
dUpp/dT 

-1.0  x  10-3  [3] 

Negative  electrode  temperature  coefficient  (VK-1), 
dUn,Q/dT 

-6.3  x  10”4  [3] 

Table  2 


Average  heat  conductivities  of  battery  components 


Negative  electrode 

Positive  electrode  Separator 

Am  (Wm-'K-1) 

1.16a 

1.14a 

0.22b 

^(porosity) 

0.25c 

0.30c 

0.74a 

ifCWm'1  KA1) 

0.57d 

0.57d 

0.57d 

&i  (fraction) 

0.29 

0.45 

0.27 

(Wn-r1  KT1) 

1.01e 

0.97e 

0.48e 

a  Ref.  [9]. 
b  Ref.  [14]. 
c  Measured  results. 
d  Ref.  [15]. 

e  The  average  thermal  conductivities  of  each  component  are  estimated  by 
Ref.  [15]. 

The  parameters  used  in  the  thermal  modeling  of  the  cylin¬ 
drical  Ni/MH  battery  are  listed  in  Tables  1  and  2. 


3.  Experimental  details 

Cylindrical  8-Ah  Ni/MH  batteries  fabricated  in  our  labora¬ 
tory  were  adopted  to  test  the  temperature  distribution  during 
charge  and  discharge.  An  8-Ah  cylindrical  Ni/MH  battery 
was  constructed  by  winding  a  positive  electrode  and  a  nega¬ 
tive  electrode  and  a  porous  polypropylene  separator  switched 
in  between  the  electrodes.  The  case  encasing  the  electrode 
groups  is  made  of  stainless  steel  with  a  thickness  of  0.5  mm. 
|3-Ni(OH)2  is  active  material  of  positive  electrode,  AB5  alloy 
powder  as  active  material  of  negative  electrode,  consists  of 
MmNi3.75Coo.55  Alo.3Mno.4.  KOH  solution  (30  wt%)  acted  as 
electrolyte  in  this  work.  The  battery  parameters  are  listed  in 
Table  1. 

Arbin  BT  Cycler  was  used  to  control  the  process  of  charg¬ 
ing  under  the  constant  current  condition.  Thermocouples 
were  used  to  measure  the  temperature  changes  on  the  bat¬ 
tery  surface.  The  data  of  current,  voltage  and  temperature 
were  recorded  by  Arbin  BT  Cycler. 

In  order  to  estimate  the  unknown  heat  capacity,  smaller 
size  AA  1.8  Ah  cylindrical  Ni/MH  batteries  with  the  same 
chemistry  and  design  were  measured  to  test  the  tempera¬ 
ture  and  voltage  curves  during  the  relaxation  period.  The 
advantage  choosing  small  batteries  is  that  their  Biot  number 
[12]  can  be  so  small  that  there  are  no  significant  tempera¬ 
ture  profiles  inside  the  batteries  under  the  normal  operating 
conditions. 


4.  Results  and  discussion 

4.1.  Estimation  of  heat  capacity 

To  estimate  the  actual  heat  capacity,  a  small  cylindrical  AA 
1.8- Ah  Ni/MH  battery  with  the  same  chemistry  and  design 
was  charged  at  1C,  2 C,  4C  rates  separately.  Since  there  is 
no  significant  temperature  profile  inside  the  small  battery 
because  of  the  small  Biot  number  [12]  under  the  condition 
given  by  Eq.  (10),  the  assumption  of  a  uniform  temperature 
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profile  across  the  small  battery  is  valid  and  hence  a  lump- 
parameter  model  of  energy  balance  can  be  applied.  The  heat 
energy  equation  (Eq.  (1))  is  thus  reduced  to  Eq.  (11), 

hL 

Bi  =  —  «  1  (10) 

A 

where  h  is  the  free-convective  heat-transfer  coefficient,  L  rep¬ 
resents  the  thickness  of  battery,  respectively 

d(CpT) 

M^-f^  =  Ah(T0O-T)+Q  (11) 

dr 

where  M  denotes  the  mass  of  battery,  A  the  surface  area 
through  which  heat  is  removed  from  the  battery  and  Q  is 
the  heat  generation  rate  in  watts. 

Cp  is  a  function  of  temperature  by  fitting  the 
temperature-time  curves  during  the  relaxation  periods,  and 
Cp  is  treated  to  be  varied  linearly  with  temperature  as  given 
by 


according  to  the  following  equations  [5], 

1  ^ — >  S,  <5P  <5S  <5n 

series-connected  —  =  >  —  =  —  H — LH -  (14) 

V 

Parallel-connected  Xz  =  =  <5p/fp  +  SsXs  +  SnXn 

(15) 

Where  <5p,  <5S  and  Sn  denote  the  volume  fractions  of  the 
positive  electrode,  the  separator  and  the  negative  electrode, 
respectively;  A.p,  Xs  and  Xn  represent  the  average  heat  con¬ 
ductivities  of  positive  electrode,  separator  and  negative  elec¬ 
trode,  respectively.  Table  2  summarizes  the  average  heat 
conductivities  of  battery  components.  According  to  Table  2, 
Eqs.  (14)  and  (15),  Xr  and  Xz  are  estimated  to  be  0.74  and 
0.85  W  m-1  K-1,  respectively,  Ref.  [9], 

Xj  =  Xm(  1  ~  s)  +  XfE  (16) 


Cp  =  Cp,ref[l  +  m(T  -  rref)]  (12) 

Where  CPiref  is  the  reference  heat  capacity  when  the  temper¬ 
ature  equals  7’ref. 

Substituting  Eq.  (12)  into  Eq.  (11),  Eq.  (13)  can  then  be 
obtained  when  the  current  is  turned  off  after  charging  finished 
[13],  i.e.  <2  =  0. 

— — — t  =  (1  +  m Too)  In  (  T°~Too\  +  2 m(T0  -  T) 
MCpM  \T-Too) 

(13) 

Cp  and  m  can  be  obtained  by  fitting  the  temperature-time 
curves  as  Eq.  (13)  measured  during  the  relaxation  periods 
following  charge.  From  Fig.  2,  Cp  ref  and  m  values  were 
obtained,  Cp  ref  =  1882  J  kg-1  K-1  and  m  =  —0.001. 


Where  s  is  porosity,  /,nl  and  Xp  are  the  heat  conductivities  of 
structural  materials  and  filling  electrolyte  materials,  respec¬ 
tively. 


4.3.  The  average  surface  temperature  profiles 


FEMLAB®,  a  commercial  FEA  software  system,  is 
employed  to  solve  the  heat  equations  in  this  work.  Fig.  3 
shows  the  charging  curves  of  8- Ah  Ni/MH  battery  at  1C,  2C 
and  4C  rates,  respectively.  The  voltage  curves  are  then  fitted 
to  calculate  the  joule  heat  in  Eq.  (1).  Open  circuit  voltage 
(OCV)  can  be  written  as  a  function  of  state-of-charge  (0) 
denoting  the  average  battery  state-of-charge  defined  to  be 


it 


Cbattery 


(17) 


4.2.  Estimation  of  heat  conductivities 

The  heat  conductivities  Xr  and  Xz  as  constants  are  highly 
anisotropic  and  microstructure-dependent  and  are  estimated 


where  t  is  the  time  in  hours,  and  Cbattery  is  the  battery  nominal 
capacity. 

The  calculated  and  measured  average  surface  temperature 
profiles  were  compared  in  Fig.  4  at  varying  charging  times 


Fig.  2.  Estimation  of  the  Cp  and  m  by  fitting  the  measured  values  during  the 
relaxation  periods. 


Fig.  3.  Charging  curves  of  an  8-Ah  cylindrical  Ni/MH  battery  during  1C, 
2 C  and  AC  constant  current  charge. 


596 


J.  Shi  et  al.  /  Journal  of  Power  Sources  157  (2006)  592—599 


Fig.  4.  The  measured  and  calculated  average  surface  temperature  variations 
of  the  battery  during  charge  at  1C,  2 C  and  4C  constant  current  under  forced 
convection  (h  =  25  W  m-2  K). 

under  the  forced  convection  (/;  =  25  W  m~2  K_1).  The  steeper 
parts  of  three  cures  mean  overcharge.  The  average  surface 
temperature  rises  from  24.4  to  29,  33  and  41  °C  at  the  end  of 
1 C,  2 C  and  4 C  charge,  respectively,  when  the  imputed  charge 
equals  110%  of  nominal  capacity,  the  average  surface  tem¬ 
perature  rises  rapidly  to  37,  42  and  51  °C  at  1C,  2C  and  4C 
overcharge,  respectively.  It  is  showed  that  the  average  surface 
temperature  increases  when  the  charging  rates  increase.  The 
battery  temperature  increase  with  time  under  all  cases,  indi¬ 
cating  that  the  battery  is  exothermic  during  charge.  When  the 
battery  was  charged  at  a  rate  as  large  as  4C,  the  calculated 
temperature  is  slightly  lower  than  the  measured  values  before 
100%  capacity.  This  is  probably  own  to  a  low  estimation  of 
the  oxygen  recombinant  heat  which  is  highly  larger  than  the 
ohmic  heat,  because  the  oxygen  reaction  current  increases 
and  the  primary  reaction  current  decreases  when  charged  at 
higher  rates.  The  results  from  the  simulation  are  well  in  agree¬ 
ment  with  experimental  measurements  at  the  range  from  1 C, 
2 C  to  AC  charge  in  spite  of  a  slight  higher  predicted  value 
during  the  overcharge.  The  rate  of  heat  generation  was  over¬ 
estimated  during  overcharge  in  Eq.  (4).  It  is  because  that  the 
oxygen  generated  at  the  positive  electrode  was  not  entirely 
reduced  at  the  negative  electrode,  but  only  a  portion  of  oxygen 
was  reduced  under  a  higher  battery  internal  pressure,  there¬ 
fore,  the  oxygen  recombinant  heat  is  not  so  much  as  that  in 
Eq.  (4). 

4.4.  Effect  of  heat-transfer  coefficients  on  temperature 
distribution 

Fig.  5  compares  the  predicted  temperature  distributions 
in  the  centre  of  the  battery  under  free,  forced  convection 
and  adiabatic  conditions.  Five  thermal  conditions  were  con¬ 
sidered  in  terms  of  an  equivalent  heat-transfer  coefficient 
imposed  on  the  surface.  Where  h  =  0  and  li  =  oo  correspond 
to  adiabatic  and  isothermal  conditions,  respectively,  h  =  6 
and  25  W  m-2  K- 1  refer  to  typical  values  of  air-free  convec¬ 


Fig.  5.  Calculated  temperature  rise  in  the  centre  of  a  battery  with  different 
heat-transfer  coefficients  (h  =  0,  6,  25,  100,  oo,  W  m~2  K_1). 

tion  and  forced  convection,  /;  =  100  Wm-2  K_1  represents 
the  extended  forced  convection. 

It  is  indicated  that  the  battery  temperature  is  remarkably 
lower  under  forced  convection  than  free  convection.  This  dis¬ 
crepancy  is  ascribed  to  the  fact  that  forced  convection  can 
sufficiently  dissipate  the  generated  heat.  The  adiabatic  behav¬ 
ior  was  also  investigated  with  the  assumption  of  the  same 
charging  characteristics;  temperature  rise  under  the  adiabatic 
condition  is  significantly  higher  than  those  under  the  con¬ 
vection  condition.  Adiabatic  modeling  is  important  because 
we  must  consider  the  consequences  of  no  heat  removal,  for 
example,  the  accident  if  the  temperature-control  system  fails. 
Under  the  adiabatic  condition,  the  temperature  increases  up 
to  about  47  °C  at  100%  charge  input,  and  dramatically  up  to 
about  60  °C  at  1 10%  charge  input,  it  is  obvious  that  adiabatic 
charge  is  highly  dangerous.  As  the  heat-transfer  coefficients 
increase,  the  temperature  in  the  centre  of  battery  is  prone  to 
a  lower  value,  when  the  heat-transfer  coefficients  increase 
close  to  infinite,  which  means  the  temperature  on  the  surface 
of  case  wall  was  constantly  kept  close  to  the  initial  temper¬ 
ature,  the  temperature  in  the  centre  of  battery  is  about  10  °C 
lower  when  compared  with  the  adiabatic  condition.  A  larger 
heat-transfer  coefficient  corresponds  to  a  larger  rate  of  heat 
dissipation,  resulting  in  a  smaller  temperature  rise  during  bat¬ 
tery  charging.  When  the  heat  transfer  coefficient  is  equal  to 
6  W  m-2  K“ 1 ,  the  battery  centre  temperature  rises  up  to  58  °C 
at  1 10%  of  the  nominal  battery  capacity.  As  the  heat-transfer 
coefficient  increases  (25  W  irT 2  K-1),  the  final  battery  tem¬ 
perature  decreases.  The  calculated  results  indicate  increasing 
heat  transfer  coefficient  can  decrease  the  centre  temperature 
which  is  the  maximum  temperature  inside  the  whole  battery. 
A  lower  temperature  can  be  obtained  at  a  higher  heat  dissi¬ 
pation  rate. 

Fig.  6  is  a  two-dimensional  temperature  profile  in  the  r- 
and  "-directions  under  the  free  (A)  and  forced  convection 
(B)  at  the  end  of  4C  charge  (110%  of  nominal  capacity).  It 
can  be  seen  that  the  average  temperature  of  the  battery  is 
strongly  affected  by  cooling  conditions.  The  higher  the  heat 
dissipation  rate,  the  lower  the  battery  temperature.  It  also 
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Fig.  6.  Temperature  profiles  in  r-  and  z-direction  at  the  end  of  4 C  charge  (A: 
/z  =  6;  B:  /*  =  25,  To  =  25  °C). 

reveals  that  the  temperature  distribution  at  the  end  of  charge 
become  more  non-uniformed  across  the  battery  especially  in 
/■-direction,  because  the  heat  is  transferred  faster  at  the  side 
of  battery  compared  to  that  in  the  centre  of  battery  when 
heat-transfer  coefficient  increases.  Although  the  temperature 
decreases  remarkably  with  a  higher  heat-transfer  coefficient, 
the  temperature  profile  across  the  battery  become  more  non- 
uniformed,  it  is  clearer  in  Fig.  7.  Obviously,  in  order  to  prevent 
the  significant  increase  of  the  temperature  in  the  battery  under 
higher  charge  rates,  the  improvement  in  cooling  conditions  is 
necessary;  however,  it  is  not  helpful  to  uniformity  of  temper¬ 
ature  profile  across  the  battery.  Instead,  it  may  cause  steeper 
temperature  profiles.  The  large  temperature  gradients  in  re¬ 
direction  may  suggest  that  the  low  conductivity  in  /‘-direction 
will  limit  heat  removal  rate,  and  therefore  the  heat  transport 
behavior  in  the  battery  is  mainly  characterized  by  the  low 
thermal  conductivity  of  separator. 

Fig.  7  depicted  the  effect  of  different  heat-transfer  coef¬ 
ficient  on  the  difference  between  maximum  and  minimum 
temperature  obtained  from  the  thermal  model  for  4 C  charge. 


Fig.  7.  Effect  of  heat-transfer  coefficients  on  maximum  and  minimum  tem¬ 
perature  under  4 C  charge. 

When  approaching  the  overcharge  region,  the  difference 
between  maximum  and  minimum  temperature  increases  dra¬ 
matically  with  a  higher  heat-transfer  coefficient.  When  the 
heat-transfer  coefficient  equals  6  W  m-2  K-1,  the  difference 
between  maximum  and  minimum  temperature  across  the  bat¬ 
tery  is  only  about  3°C  at  the  end  of  4C  charge,  while  the 
heat  transfer  coefficient  increases  up  to  100  Wm”2  K-1,  the 
value  is  remarkably  up  to  about  15  °C  at  the  same  condi¬ 
tion.  Non-uniformed  temperature  profiles  result  in  different 
electrochemical  behaviors,  and  in  turn,  cause  the  battery  per¬ 
formance  fade.  Together  with  Fig.  6,  it  is  shown  that  the 
only  increase  of  the  heat  transfer  coefficient  cannot  efficiently 
improve  the  uniformity  of  temperature  across  the  whole  bat¬ 
tery  in  spite  of  a  lower  average  temperature  the  system  can 
obtain.  It  can  also  be  seen  that  the  worse  results  encounter 
when  overcharged,  the  temperature  differences  trend  to  larger 
values  when  the  heat-transfer  coefficients  increase. 

4.5.  Effect  of  heat  conductivities  on  temperature  profile 

The  heat  conductivity  is  decided  by  the  battery  struc¬ 
ture,  electrode  structures  and  materials.  According  to  Ref. 
[9],  heat  conductivity  is  related  to  the  electrode  porosity  and 
the  ratio  of  positive,  negative  electrode  and  separator,  which 
decide  the  electrochemistry  of  battery  system;  therefore,  it  is 
not  practical  to  change  the  heat  conductivity  in  large  scale 
for  the  actual  battery.  To  study  the  effect  of  heat  conduc¬ 
tivities  on  battery  thermal  behavior,  an  assumption  with  a 
small  change  range  of  heat  conductivity  was  adopted  in  this 
work.  Fig.  8  compared  the  effect  of  different  heat  conductiv¬ 
ity  on  the  battery  temperature  profile  in  /  -direction.  It  can  be 
seen  that  the  maximum  temperature  decreases  and  minimum 
temperature  increases  with  the  increase  in  thermal  conduc¬ 
tivity  in  /--direction.  When  Xr  rise  about  30%  from  0.74  to 
0.96  W  m_  1  K_  1 ,  although  the  battery  temperature  decreases 
slightly,  the  temperature  profile  becomes  more  uniform,  the 
difference  between  the  maximum  and  minimum  temperature 
across  the  battery  decreases  from  3.5  to  2.7  °C,  while  drops 
about  30%  from  0.74  to  0.56  W  m_  1  K- 1 ,  the  battery  temper¬ 
ature  increases  slightly,  but  the  temperature  profiles  become 
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Fig.  8.  Temperature  profiles  under  different  conductivities  in  /--direction 
(z=Ht 2)  at  the  charge  end  of  4 C  rate  (Xz  =  0.84  Wm  1  K-1,  Tq  =  25  °C). 

more  non-uniformed,  the  difference  between  the  maximum 
and  minimum  temperature  across  the  battery  rises  from  3.5 
to  4.6  °C.  It  seems  that  the  increase  of  thermal  conductivity 
can  obtain  a  more  uniformed  temperature  profile  based  on  the 
thermal  modeling.  In  general,  the  increase  of  the  heat  con¬ 
ductivity  can  improve  the  uniformity  of  battery  temperature 
profile  to  some  extend  in  the  actual  battery,  but  this  leads 
probably  to  the  battery  electrochemistry  changes  which  con¬ 
tradict  to  the  improvement  of  heat  conductivities.  However, 
from  a  heat  transport  point  of  view,  some  better  designs  could 
be  achieved  by  which  temperature  gradient  within  a  battery 
is  minimized  [5]. 

4.6.  Effect  of  charging  rates  on  temperature  profiles 

Fig.  9  shows  temperature  profiles  of  an  8-Ah  Ni/MH  bat¬ 
tery  in  /-direction  (z  =  H/ 2)  under  different  charge  rates.  It 
can  be  seen  that  under  low  charge  rate,  the  temperature  in 
the  battery  presents  a  small  increase,  however,  under  high 
charge  rates,  the  temperature  rise  in  the  battery  is  very  high. 
The  maximum  temperature  in  the  battery  rises  from  25  to 
31,  36  and  45  °C  at  the  end  of  1C,  2C  and  4C  charge, 
respectively.  Furthermore,  when  charge  rates  increase,  the 
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Fig.  9.  Temperature  profile  in  the  battery  ( z  =  HI 2)  in  r-direction  under  differ- 
ent  charge  rates  (\r  =  0.74 Wm-1  K“\  =  0.84 WirT1  K-1,  7b  =  25°C, 
/i  =  25  W  m-2  K-1). 

temperature  gradients  of  the  battery  become  steeper  under 
the  same  condition.  The  difference  between  maximum  and 
minimum  temperature  rises  from  0  to  3,  5  and  7  °C  at  the  end 
of  1 C,  2C  and  4C  charge,  respectively.  The  similar  results 
can  be  obtained  in  "-direction.  The  calculated  results  indi¬ 
cate  that  a  higher  charge  rate  leads  to  a  higher  temperature  rise 
and  a  less  uniformed  temperature  profile,  which  are  harmful 
to  the  battery  performance.  Obviously,  a  more  efficient  heat 
removal  should  be  selected  in  order  to  adopt  higher  charge 
rates,  which  results  in  much  faster  heat  generation  rate. 

4.7.  Effect  of  overcharge  on  temperature  profiles 

Fig.  10  shows  the  effect  of  overcharge  on  the  tempera¬ 
ture  profiles  of  the  battery.  It  can  be  seen  that  the  overcharge 
results  in  a  serious  temperature  rise  and  a  less  uniformed  tem¬ 
perature  profile.  Almost  all  the  current  applied  to  the  battery 
is  used  to  generate  oxygen  at  the  positive  electrode  during 
overcharge.  The  generated  oxygen  at  the  positive  electrode 
can  be  reduced  at  the  negative  electrode,  and  this  reduction 
reaction  will  generate  much  more  heat  than  that  from  bat¬ 
tery  internal  resistance,  and  hence  a  steeper  temperature  rise. 


(B)  r(m) 


Fig.  10.  Effect  of  overcharge  on  the  temperature  profiles  in  r-direction  at  the  charge  end  of  2 C  rate  (A)  and  4 C  rate  (B). 
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When  input  charge  rises  10%  from  90%  to  100%  of  nomi¬ 
nal  capacity,  the  maximum  temperature  rises  only  1  and  2  °C 
under  2 C  and  4 C  charge,  respectively,  while  input  charge  rises 
10%  from  100%  to  1 10%  of  nominal  capacity,  the  maximum 
temperature  increasingly  rises  10  and  12  °C  under  2 C  and 
4 C  rates  overcharge,  respectively.  The  difference  between 
maximum  and  minimum  temperature  is  5  and  7  °C  at  100% 
charge  input  under  2 C  and  4 C  rate  charge,  respectively,  while 
the  difference  rises  even  up  to  7  and  9  °C  at  10%  overcharge 
under  2 C  and  4 C  rates,  respectively.  Obviously,  overcharge 
leads  to  overheat  and  hence  a  stronger  temperature  rise  and 
a  steeper  temperature  gradient.  Therefore,  overheat  can  lead 
to  a  violent  change  of  battery  electrochemistry,  which  will 
result  in  dangerous  situation.  Overcharge  modeling  can  give 
some  important  information  of  system  temperature  variations 
when  charge-control  system  fails.  The  calculation  indicates 
that  it  is  necessary  to  avoid  overcharge  when  fast  charging 
the  Ni/MH  battery. 

5.  Conclusions 

Thermal  analysis  of  cylindrical  nickel-metal  hydride 
batteries  during  fast  charge  has  been  carried  out  by  the 
two-dimensional  thermal  modeling  proposed  in  this  work. 
Detailed  information,  such  as  effects  of  dissipating  heat 
coefficient,  heat  conductivity,  charge  rate  and  overcharge 
on  the  battery  temperature  profile  across  the  battery  during 
fast  charge  can  be  digitally  obtained.  The  thermal  behav¬ 
ior  of  batteries  can  be  better  understood  with  these  simulated 
results  from  the  mathematical  modeling.  Based  on  the  results 
obtained  from  model  calculation  and  experimental  measure¬ 
ments,  it  can  be  seen  that  the  increase  of  the  heat-dissipating 
coefficient  can  decrease  the  battery  temperature,  but  leading 
to  a  less  uniformed  temperature  profile  across  the  whole  bat¬ 
tery.  A  more  uniformed  temperature  profile  can  be  obtained 
when  adopting  a  large  heat  conductivity,  which  may  probably 
be  achieved  by  improving  the  battery  electrode  structure  and 


materials.  When  charge  rates  increase  the  battery  temperature 
increasingly  rises  and  the  temperature  profiles  become  less 
uniformed  under  the  same  input  charge,  overcharge  results 
in  a  dangerous  temperature  rise,  which  should  be  avoided 
carefully. 
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